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' The element Vub of the CKM mixing matrix is the smallest of the quark cou- 

CN ! plings and is crucial in understanding CP violation in the B system. In this review, 

' we discuss the present status of the determination of its magnitude, which involves 

. a significant effort both experimentally and theoretically. Decisive progress has 

I been achieved in recent years, thanks in part to the large data sets that have be- 

p • come available at the B factories. Based on these data we propose an average |Kib| 

^ , value of (3.67 it 0.47) x 10"'^, which has an uncertainty of 13 %. With the experi- 

Qj I mental and theoretical progress expected in the next few years a determination to 

• an accuracy of 10 %, or better, seems feasible. 
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To appear in Review of Particle Properties 

The precise determination of the magnitude of V^b with a robust, well-understood uncertainty 
remains one of the key goals of the heavy flavor physics programs, both experimentally and the- 
oretically. Because \Vub\, the smallest element in the CKM mixing matrix, provides a bound on 
the upper vertex of one of the triangles representing the unitarity property of the CKM matrix, it 
plays a crucial role in the examination of the unitarity constraints and the.^fundamental questions 
on which the constraints can bear (see the minireviews on the CKM matrixtl and on CP-violationl3 
for details). Investigation of these issues require measurements that are precise and that have 
well-understood uncertainties. 

The charmless semi-leptonic (s.l.) decay channel b ulu provides the cleanest path for the 
determination of |Vuf,|. However, the theory for the heavy-to-light b ^ u transition cannot be 
as well constrained as that fox the heavy-to-heavy b ^ c transition used in the determination of 
\Vcb\ (see the \Vcb\ minireviewB) . The extraction of \Vub\ and the interplay between experimental 
measurements and their theoretical interpretation are further complicated by the large background 
from b — > (dv decay, which has a rate about 60 times higher than that for charmless s.l. decay. 
Measurements based both on exclusive decay channels and on inclusive techniques have been, and 
are being, pursued. 
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The last several years have seen significant developments in both the theoretical framework 
and the experimental techniques used in the study of 6 — > uii'. The inclusive theory has pro- 
gressed significantly in categorization of the corrections to the base theory still needed, in their 
relative importance in different regions of phase space, and in the determination of some of them. 
Recent work on exclusive processes bolsters confidence in the current uncertainties for the form 
factor calculations needed to extract \Vub\- Experimentally, we have new inclusive and exclusive 
measurements that minimize dependence on detailed modeling of the signal process to separate 
signal from the b civ background, have well-defined sensitivities in particular regions of phase 
space and have improved signal-to-background ratios. These improvements provide us with a first 
opportunity to develop a method for a robust determination of with complete error estimates, 
including constraints on hitherto unquantified contributions. We review the current determinations 
of \yub\i focusing primarily on these recent developments. An average of the inclusive information 
from all regions of phase space remains, unfortunately, beyond our reach because of the potentially 
sizable corrections for which we lack estimates. Rather, we combine the inclusive results to obtain 
a central value and, in particular, a more complete evaluation of the uncertainty than has been 
possible in the past. 



Inclusive measurements of 6 — ulv 



Theoretically, issues regarding the calculation of the total semilentonic partial width Y{B ^ X^iv^ 
via the operator product expansion (OPE) are well-understoodaEl. The OPE is both a nonper- 
turbative power series in l/mt and a perturbative expansion in a^. At order it predicts 
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where A2 parameterizes the hyperfine interaction between the heavy quark and the light degrees of 
freedom and Ai is related to the Fermi momentum of the heavy quark. The perturbative corrections 
are known to order agEHl. The OPE is alternatively expressed in terms of the nonperturbative 
parameters //^ and //g, which are closely related to Ai and A2, respectively, but differ significantly 
in their infrared treatment. Within the OPE, the importance of a proper field theoretic treatment 
of the parameters is paramount both for the total rate and for the restricted phase space studies 
discussed below. The treatment of the quark mass and its associated uncertainties are particularly 
important given the strong mass dependance of the width. Such considerations have led to useful 
definitions like the kinematic mass, which are discussed in detail in Ref. 11 and Ref. 12| . 

The error indiLced by uncertainties in the nonperturbative parameters Ai_2 is relatively small, 
and an evaluationlij by the LEP VUB working group yielded 
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The quoted uncertainty is dominated by the uncertainty in the b quark mass, for which m\f {IGeV) = 
4.58 ibJ1.09 GeV was assumed. The value and uncertainty are in good agreement with a recent 
surveytj. No weak annihilation uncertainties (discussed below) are included in the quoted OPE 
error. Use of the quark-level OPE for prediction of moments of the true inclusive spectra has gen- 
erated concern regarding potential violation of the underlying assumption of global quark-hadron 
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duality. This concern has been confronted by theoretical wisdomEl supporting global duality for 
the inclusive b uli' transition, and new data both support this assumption and allow placement 
of quantitative limits on the violation. In particular, the exclusive and inclusive extractions of \Vcb\ 
agree to (0.8 ± 1.6) x lO^^B. Taking the uncertainty as an upper bound on the global duality 
violation for \Vub\ shows that it should not exceed ~4%. Using this as bound as an uncertainty 
estimate for duality effects in the partial width prediction brings the total uncertainty on \Vub\ to 
6.8%. The agreement of the OPE parameters extracted using moments of different distributions in 
s.l. decays further supports the small scale for duality violation effects. 

While theoretically the total inclusive rate would allow determination of \Vub\ to better than 
10%, experimentally the much more copious b ciu process makes a measurement over the full 
phase space unrealizable. To overcome this background, inclusive b uiu measurements utilize 
restricted regions of phase space in which the b civ process is kinematically highly suppressed. 
The background is forbidden in the regions of large charged lepton energy > {Mg — M^) / {2Mb) 
(the endpoint), low hadronic mass Mx < Mu and large dilepton mass > {MB — M^y. Extraction 
of \Vub\ from such a measurement requires knowledge of the fraction of the total b — > uiv rate 
that lies within the utilized region of phase space, which complicates the theoretical issues and 
uncertaintv-considerahly. Ref. 16 and |Ref. 17 discuss the issues in detail. 

CLEOtS, BaBarcJ and BELLEE3 have all presented recent measurements of the b u£u rate 
near the endpoint. The results, which are for integrated ranges in the T{4S) rest frame, are sum- 



marized in p?able 1.1| . Experimentally, these measurements must contend with a large background 
from continuum e+e~ annihilation processes. Suppression of these backgrounds introduces signif- 
icant efficiency variation with the of the decay, which introduces model dependence. Greater 
awareness of this issue has resulted in more sophisticated suppression methods in these recent mea- 
surements and thus in over a factor of three reduction in the model dependence of the measured 
rates relative to earlier measurementsEilE3. Future measurements, either using fully-reconstructed 
B tag samples that would remove the problem or a modest binning, would essentially eliminate 
the remaining model dependence. 



Table 1.1: Partial branching fractions for b uiu within the charged lepton 
momentum range (T(45') frame) from 2.6 GeV/c down to the indicated minimum. 
The estimated fraction /e of the total b — uii' rate expected to lie in that range is 
also given. The dagger (|) indicates the quantity that received the QED radiative 
correction appropriate to the indicated mix of electrons and muons, which has not 
always been treated self-consistently in the literature. 

(GeV/c) (10^^) 



2.0 4.22 ±0.33 ±1.78 tO.266 ± 0.041 ± 0.024 CLEO(e,/i) 

2.1 3.28 ± 0.23 ± 0.73 t0.198 ± 0.035 ± 0.020 CLEO(e,/i) 

2.2 2.30 ±0.15 ±0.35 t0.130 ± 0.024 ± 0.015 CLEO(e,/i) 

2.3 1.43 ±0.10 ±0.13 t0.074 ± 0.014 ± 0.009 CLEO(e,/i) 
tl.52 ± 0.14 ± 0.14 0.078 ± 0.015 ± 0.009 BaBar (e) 
1.19 ±0.11 ±0.10 t0.072 ± 0.014 ± 0.008 BELLE (e) 

2.4 0.64 ±0.07 ±0.05 t0.037 ± 0.007 ± 0.003 CLEO(e,/i) 



BaBaiB and BELLEEd have presented new analyses of the low Mx regionc^TEa. They also 
utilize a moderate (only ~ 10% loss) pi > 1.0 GeV/c requirement. This technique was pioneered by 
Delphi at LEPE3, but the achievable resolution and signal-to-background ratio were lower compared 
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to those obtained at the B factories. Because of experimental resolution on Mx, the b ciu 
background smears below its theoretical lower limit of Mx = M^, so experiments must impose 
more stringent Mx requirements which are theoretically more problematic. The BaBar and BELLE 
analyses are based on "i? tag" samples of fully reconstructed hadronic decays and D^*Hv decays, 
respectively. Li both cases, Mx is calculated directly from the particles remaining after removal 
of tag and lepton contributions. The BaBar analysis, in particular, reveals a beautiful h uii' 
signal with an unsurpassed signal to background ratio of about 2:1 in the region Mx < 1.55 GeV/c, 
which rivals that of current exclusive analyses. This analysis demonstrates the anticipated power 
of a large fully reconstructed sample, both in the signal to background levels and in the excellent 
resolution that can be achieved. The efficiency versus Mx appears reasonably uniform, and the 
signal yield fitting procedure minimizes the dependence of the extracted rate on the modeling of 
the detailed b uiv dynamics. Both allow for improved determination of \Vub\ as theory advances. 

Determination of the fraction of the b uiu rate in the a^_.eiidpoint or the low Mx region 
requires resummation of the OPE to all orders in Ex^qcd/itT'x^^^"- The resummation results, 
at leading-twist order, in a nonperturbative shape function f{k^), where = k'^ + k\\ and k^ = 
p'^ — nihV^ is the residual b quark momentum after the "mechanical" portion of momentum is 
subtracted off. Spatial components k\\ and k± are defined relative to the nibv'^ — (roughly the 
recoiling u quark) direction. At this order, effects such as the "jiggling" of kj^ are ignored and 
the differential partial width is given by the convolution of the shape function with the parton 
level differential distribution. Because the shape function depends only on parameters of the B 
meson, this leading order description holds for any B decay to a light quark. It holds, in particular, 
for B ^s£i^ which can provides an estimation of f{k^) via the shape of the photon energy (E^) 
spectrumEilO. In addition to the increased uncertainty on \Vub\ from the m?, and b quark kinetic 
energy contributions that results from the restridupn of phase space, higher twist contributions and 



unknown power corrections of order Aqcd/-^sSS also contribute to the uncertainty, as will be 
discussed further below. 

Ideally, \Vub\ would be determined without introductianof an intermediate extracted shape 
function through the use of appropriately weighted spectraHl£ZrO. This would avoid introduction 
of an element of model dependence. For the lepton spectrum, for example, one would take 



V, 



ub 



VtbV,* 



- ^pert£4#T + O (AqCd/M^) , (1.3) 



where Kp^rt is a calculable perturbative kernel, abd r„(£'o) and Ts{Eo) are appropriately weighted 
integrals over, respectively, the lepton energy and photon energy spectra above the minimum cutoff 
energy Eq. Practical application of this approach awaits measurement of the lepton momentum 
spectrum in the B, not the T(4S'), rest frame, which B tag methods willjj)£rmit in the future. 
Similar expressions exist for integration over the low hadronic mass regioncScJ, so, in principle, 
current Mx analyses should be able to take such an approach. Experimental efficiency and lepton 
momentum cutoffs must, however, be incorporated into the integrals. To date, experiments have 
instead introduced intermediate shape functions of a variety of forms into the inclusive analyses, 
as discussed below. 

A third way to isolate the charmless s.l. signal is to use a selection based on the q'^ of the 
leptQnici_system. Restriction of phase space to regions of large q^ also restores the validity of the 
OPEBH and suppresses shape function effects. Taking only the region kineinatically forbidden 
to 6 ^ ciu, q^ > {Mb — M^Y, unfortunately introduces a low mass scaleEjIla into the OPE and 
the 1/m^ uncertainties blow up to be of order (Aqcd/'it^c)^ ■ However, a combination of Mx with 
looser q^ requirements can suppress both b — > ciu background experimentally and shape function 
effects theoretically. Furthermore, the q^ requirement moves the parton level pole away from the 
experimentally feasible Mx requirement. The shape function effects, while suppressed, cannot be 
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neglected. One drawback of the requirement is the elimination of higher energy hadronic final 
states, which may exacerbate duality concerns. 

A recent BELLE analysii3 has been performed in this region. BELLE employs a > 1.2 
GeV/c requirement in the T(45') rest frame, and an "annealing" procedure to separate recon- 
structed particles into signal and "other B" halves. They then examine the integrated rate in the 
region Mx < 1-7 GeV and > 8 GeV^ to extract |Kt;,|, which again has the desired effect of 
minimizing dependence of the analysis on detailed h — > ulv modeling The signal to background 
ratio of the annealing technique, about 1:6 for the BELLE analysis, is significantly degraded rela- 
tive to that of the hadronic B tag technique. As we mentioned in the previous review, control of 
background subtractions of this size requires extreme care and careful scrutiny of the associated 
systematic issues. BELLE finds the rate Ai3 in that region of phase space to be 

AB = (7.37 ± 0.89,tat ± 1.12sy. ± 0.55,,, ± 0.24„,,) x 10"^ 

An analysis of this jcestricted region of phase space, for which the shape function influence is 
signiflcantly reducedcS, with the signiflcantly cleaner B tag technique should be a priority for both 
B experiments. 

Each analysis discussed here has relied on an intermediate shape function to evaluate the 
fraction of the inclusive rate that lies in its restricted region of phase space. The endpoint analyses 



have used rate fractions Ref. 18 based on intermediate shane functions derived from the CLEO 
b photon spectrum. Several two-parameter ansaetzeElrE3, F[A'^^, Af''^], were implemented 

as the form of the shape function. These parameters are related to the HQET parameters of 
similar name, and play a similar role in evaluation of the rates. At this time, however, we do 
not know the precise relationship between the shape function parameters, or the moments of the 
shape function, and the HQET nonperturbative parameters A and AiOt3. The fact that A'^''^ and 
Af^ depend on the functional ansatz while the HQET parameters depend on the renormalization 
scheme underscores the current ambiguity. 

With the limited statistics available, there exists a strong correlation between the two 
parameters because of the interplay between the effective b quark mass (controlled by A'^''^) and the 
effective b quark kinetic energy (controlled by Af''^) in determining the mean of the spectrum. No 
external constraints that could break the correlation, such as ml or measured b — > civ' moments, 
have been input into the b ^ s'j fits because of their unknown relationship to the shape function 
parameters. The resulting effective rrib mass range contributing to the uncertainties (±200 MeV) 
is therefore much larger than the current nii, uncertainty. Given the current independence of the 
shape function and rrib determinations and the broad effective rrih range sampled in the shape 
function, we do not consider it necessary to treat the -E^-derived phase space fractions and the 
partial width (Eq. (1.2)) as positively correlatedu^. As the data statistics increase, it will become 
possible to constrain the shape function parameters directly from distributions in 6 — > uiv', such as 
the Mx spectrum, thus removing the uncertainties introduced from their derivation from another 
class of decays. Furthermore, once the renormalization behaviour of the shape function and the 
relationship of its moments to the HQET parameters becomes known, the powerful constraints 
from the kinematic mass of the b quark and from moments information in the b — > clu system can 
be incorporated into a shape function derivation based either on 6 — > 57 or 6 ^ uii/. 

Two alternate approaches to the shane function evaluation have been taken in experimental 
studies so far. In their low-Mx analysis^, BaBar has evaluated the phase space fraction using 
the same f{k+) parameterizations noted above, but has substituted HQET parameters derived 
from stmiies of spectral moments of the 6 ^ S7 and b c^tjj)rocesses. The BELLE Mx-q^ 
analysig^( discussed below) uses the calculation of Bauer et a/.cS based on a form with a single 
parameter a = A^^ /Xf^ , which was estimated from the ml^^^ mass and from typical estimates for 
Ai. The uncertainties in the different rate fractions in the momentum endpoint, the low Mx and 
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the Mx~q^ regions are strongly correlated, and the values and uncertainties are sensitive to the 
theoretical assumptions made. Hence, a common theoretical scheme must be chosen for meaningful 
comparison of the extracted values of \Vub\- Given the ad hoc nature of the association of shape 
function parameters with the HQET parameters in the MS or the T(IS') mass schemeCl, and the 
difficulty in evaluating the uncertainty in such an association, we have chosen to extract \Vub\ from 
all of the measurements discussed using the b — > s7-derived shape function. 



Table 1.2: Summary of inclusive \Vub\ measurements. The last five measure- 
ments are incorporated into the analysis presented below. The errors in the first 
group are the experimental and theoretical uncertainties. The errors in the second 
group are from the statistical, experimental systematic, -E^-based rate fraction, 
and Ftot uncertainties. The two groups are not directly comparable as they have 
not been evaluated with identical theoretical inputs. 













\Vub\ 


:io-3) 






ALEPH [53] 








4 


12 ±0 


67 ±0 


76 


neural net 


L3 [54] 








5 


70 ± 1 


00 ± 1 


40 


cut and count 


DELPHI 








4 


07 ±0 


65 ±0 


61 


Mx 


OPAL [55] 








4 


00 ±0 


71 ±0 


71 


neural net 


LEP Avg. 








4 


09 ±0 


37 ±0 


56 




CLEO [56] 








4 


05 ±0 


61 ±0 


65 


dV/dq^dM^dEe 


BELLE 








5 


00 ±0 


64 ±0 


53 


Mx, D^*Hv tag 


CLEO 


4 


11 







13 ±0 


31 ±0 


46 ±0.28 


2.2 < p < 2.6 


BaBar 


4 


31 







20 ±0 


20 ±0 


49 ± 0.30 


2.3 < p < 2.6 


BELLE 


3 


99 


± 





17 ±0 


16 ±0 


45 ±0.27 


2.3 < p < 2.6 


BELLE 


4 


63 







28 ±0 


39 ±0 


48 ± 0.32 


Mx < 1.7, g2 > 8 


BaBar 


4 


79 







29 ±0 


28 ±0 


60 ± 0.33 


Mx < 1.55 



The full set of inclusive \Vub\ results is summaxized in p?able 1.^ , which is an updated version 
of the Heavy Flavors Averaging Group summaryH. All endpoint results have QED radiative 
corrections applied correctly. The listed uncertainties do not include contributions for potentially 
large theoretical corrections that have been categorized but remain incalculable (see below). The 
last five results in the table, which we will use below, have been updated to a common framework 
based on the CLEO ii^^-derived shape function. The rate fractional for the BaBar Mx analysis 
(/m) and the BELLE Mx - analysis /,m are fM = 0.55 ± 0.14 and f^u = 0.33 ± 0.07. The 
central values for these and for the endpoint fractions ( [Table lA ) correspond to an exponential shape 
function ansatdla and (Af^,A'^^) = (—0.342,0.545), with small corrections related to background 
subtractions in the 6 — > 57 spectrum. The errors are dominated by the statistical uncertainty 
in the f{k+) fit to the Ej spectrum, but include contributions from experimental systematics, 
uncertainties and modeling. Incorporation of results beyond those used here will require significant 
input from the experimental analyses, and is left to the HFAG. 
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1.1. Combining inclusive information 

Evaluation of the total uncertainty on \ Vub\ remains problematic because of a variety of theoretical 
complications. A recent reviewllj discusses these issues in detail. There are three main contributions. 
IDia^ first arises from subleading (higher twist) contributions to the shape function resummation 
e3~e9. These involve incorporation of effects such as the variation of k±, and are not universal for 
all B decay processes. Hence with the use of 6 — > 57 to obtain a shape function, there are two 
contributions, one from subleading contributions to the use of a shape function in 6 — > uiu process 
itself, and the second from the different corrections in 6 — > 57 from which the shape function is 
obtained. These contributions ap^ potentially large, since they are of order AQco/^b- Indeed, 
a partial estimate of these effects^ for the momentum endpoint region finds corrections that are 
similar in size to the total uncertainties of those analyses. 

The second contribution, from "weak annihilation" pmcesses, is formally of order {^qcd /"^bY 
but receives a large multiplicative enhancement of IGyr^EilS. The contribution, which requires 
factorization violation to be nonzero, is expected to be localized near ~ m^, and this localization 
can result in a further enhancement of the effect on \ Vub\- For the endpoint region, which sees about 
10% of the total rate, an effect on the total rate of 2-3% (corresponding to factorization violation 
of about 10%), produces an effect on the measured rate of 20-30%. 

Finally, there are unknown contributions from potential violation of local quark hadron duality. 
The true differential distribution cannot be predicted via the OPE ~ the resonant substructure is not 
described. However, spectra integrated over a sufficiently broad range should be better described. 

The problems just outlined present a challenge to the averaging of the various inclusive results. 
Results with a potentially large bias might be included with neither a correction nor an appropriate 
uncertainty due to these effects. The resulting \Vub\ determination would be potentially biased and 
the attached uncertainty unreliable. This motivated us not to provide an average result in the first 
edition of this review two years ago. 

As an alternative, we here choose measurements in the region of phase space that appears to 
have the best compromise of the affects discussed to to obtain an estimate of \Vub\- Measurements 
from the other regions of phase space, which have increased sensitivity to one or more of the 
corrections, then provide limits on the uncertainties from these effects and thereby allow as complete 



as possible an estimation of the theoretical uncertainty, as first proposed by Ref. 52. At this time, 
the low Mx, high region appears to be the best motivated choice. It has reduced (though 
by no means negligible) corrections from the shape function and thus also from the subleading 
contributions to the shape function. Yet it integrates over a sufficient fraction of the spectrum to 
dilute weak annihilation contributions and concerns on local quark hadron duality. 

While this choice is at present subjective, it offers the advantage of a reduction of the shape 
function influence coupled with the ability to bound the remaining theoretical uncertainties. In 
the opinion of the reviewers, this is a reasonable tradeoff for the statistical loss relative to the low 
Mx region. We expect that each experiment will perform an improved combination of information 
from the different regions of phase space where the experimental and theoretical correlations can 
be made manifest more straightforwardly. 

We further stress that we view all three regions as equally crucial in this combination of 
information, as a more complete evaluation of the inclusive uncertainty than has previously existed 
is necessary for proper use of the inclusive results. The choice of the phase space region should not 
be misconstrued as a preference of experimental technique. Indeed, we look forward to a similar 
(or improved) analysis when a sample of clean results based on fully tagged B samples have been 
obtained for all regions of phase space. 
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At present only BELLEE^ has contributed a result for this region of phase space, so for now we 
take this result as the "central value": 

I Kb 1/10-3 = 4.63 ± 0.28,tat ± 0.39,y,0.48f^,, ± 0.32rthy 

i (TwA i (TssF i (Tlqd- 

Additional measurements by the B factories of the rate in this region of phase space will soon 
improve the experimental uncertainties. 

We must determine the last three uncertainties for weak annihilation (WA), subleading shape 
function corrections (SSF) and local quark hadron duality (LQD). The measurements from other 
regions of phase space are crucial for this task. 

We assume that the WA contribution is largely contained within each of the pi > 2.2 GeV/c, 
the Mx < 1.55 GeV and the combined Mx < 1-7 GeV, > 8 GeV^ regions. The contribution 
will be most diluted in the low Mx region, with the rate fraction fm = 0.55 it 0.14, and most 
concentrated in the endpoint region, with the rate fraction /g = 0.14 it 0.03 (without radiative 
corrections). It is simple to show that for a neglected WA contribution, a comparison of \Vub\ from 
these two regions would predict the bias in the Mx, region (with rate fraction f^M = 0.33 ±0.07) 
to be 

[(1 - Ui) /Ui] [fefM/ (/m - /.)] « 0.39 (1.4) 

of the observed difference. Comparison of the endpoint result from CLEO and the low Mx result 
from BaBar, taking into consideration the almost total correlation in the shape function and Ttot 
uncertainties, yields AlT^bl/lO^'^ = 0.69=b0.53. There is not sufficient sensitivity to draw conclusions 
regarding the presence of a WA component, but we can place a bound. We take the larger of the 
error and central value and scale according to Eq. (1.4) to obtain 

0"wA ~ 0.27. 

To estimate the uncertainty from the subleading corrections to the shape function, we assume 
that subleading corrections will scale like the fractional change in the predicted rate (AT /T) with 
and without convolution of the parton-level expression with the shape function. As the base 
comparison, we take the low AIx region, with {AT/T)m = 0.15, and compare to the combined Mx, 
region, with {AT /T)qM = —0.075). The shifts again depend on the shape function modeling, and 
the quoted values correspond to the f{k^) from the best fit to the CLEO spectrum . The theory 
uncertainties are again correlated, and we find A|V^;b|/10~'^ = 0.16 ± 0.63. Scaling the uncertainty 
of the comparison by \{AT /T) qM / {^'^ l'^) m\ = 0.49, we have 

(TssF ~ 0.31. 

Finally, we must make an estimate of the local duality uncertainty. We assume that a potential 
violation will scale with the fraction of rate / in a given region as (1 — /)//. This form ranges 
from no "local" violation for integration of the full phase space (/ = 1), to large uncertainty for 
use of a very localized region of phase space (/ — > 0). The estimate derives from comparison of the 
CLEO pi > 2.2 GeV/c analysis (/ ~ 0.14 ± 0.03) to the average of the BaBar and BELLE p/, > 2.3 



GeV/c analyses (/ ~ 0.07 ± 0.02). The subleading correction estimates of |Ref. 59| are applied to 
minimize potential cancelation between duality violation and subleading corrections. This yields 
(|Kbp-3-|K6p-H0.27)/10-3 = 0.29±0.38, where the 0.27 is the estimate of the relative subleading 
correction. With our scaling assumption, we then apply a scale factor s of 

„ _ (-*■ ~ fg^') I --1 29 



(1 - /2.3) //2.3 - (1 - /2.2) //2 
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to the uncertainty in this difference estimate. Our local duality estimate therefore is 

(^LQD ~ 0.11. 

From this analysis, we finally obtain 

iKol/lO-^^ = 4.63 ± 0.28,tat ± 0.39,y, ± 0.48/,„ ± 0.32rthy 
±0.27wa±0.31ssf±0.11lqd, 

for a total theory error of 15% and total precision of 18%. Given that the uncertainties are domi- 
nated by experimental limits, addition in quadrature seems appropriate. Note that these estimates 
apply only in the combined low Mx, high region of phase space. The limits presented here can 
be improved both in robustness, through more sophisticated scaling estimates, and in magnitude, 
through additional and improved \Vub\ measurements and through inputs from other sources. The 
consistency of the values of \Vub\ extracted with different inclusive methods and the stability of 
the results over changes in the selected region of phase space will provide increasing confidence in 
the reliability of the extracted results and of their estimated uncertainties. Improvement of the 
b ^ s'j photon energy spectrum is key until a self-consistent extraction of the shape function from 
b uiv transition becomes available. Comparisons of the versus Ds semileptonic widths and 
of the rates for charged versus neutral B mesons can provide estimates of the weak annihilation 
contributionsEj. Finally, improved theoretical guidance concerning the scaling of the effects over 
phase space would allow development of a simultaneous extraction of \Vub\ and the corrections, with 
all experimental information contributing directly to \Vub\- 



2. Exclusive measurements of 6 ^ uIp 

Reconstruction of exclusive b uiv channels provides powerful kinematic constraints for sup- 
pression of the h civ background. For this suppression to be effective, an estimate of the four 
momentum of the undetected neutrino must be provided. The measurements to date have made use 
of detector hermeticity and the well-determined beam parameters to define a missing momentum 
that is used as the neutrino momentum. Signal-to-background ratios (S/B) of order two have been 
achieved in these channels. 

To extract \yub\ from an exclusive channel, the form factors for that channel must be known. 
The form factor normalization dominates the uncertainty on |V^6|. The g^-dependence of the form 
factors, which is needed to determine the experimental efficiency, also contributes to the uncertainty, 
but at a much reduced level. For example, the requirement of a stiff lepton for background reduction 
in these analyses introduces a g^-dependence to the efficiency. In the limit of a massless charged 
lepton (a reasonable limit for the electron and muon decay channels), the B — > -kIv decay depends 
on one form factor fi{q^): 

dV (B^ ^ vr-^+z.) _ , G%plMl 

dyd cose, 327r3 ^^'^^ ' ' ' 

where y = q^/M^ and 9, is the angle between the charged lepton direction in the virtual W {l + u) 
rest frame and the direction of the virtual W . For the vector meson final states p and w, three 
form factors Ai, A2 and V are necessary (see e.g. referenceE^) . 

Calculation of these form factors constitutes a considerable theoretical industxj^) -ffiith a variety 
of techniques nmv being employed. Form factors based on lattice QCD calculationd^lTLj and on light 
cone sum ruleaHTEJ currently have uncertainties in the 15% to 20% range. A variety of quark model 
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calculations existsEHTl2J. Finally, a number of other approachest2irii22, such as dispersive bounds 
and experimentally-constrained models based on Heavy Quark Symmetry, seek to improve the (f' 
range where the form factors can be estimated, without introducing significant model dependence. 

Of particular interest are the light cone sum rules (LCSR) and lattice QCD (LQCD) calcula- 
tions, which minimize modeling assumptions as they are QCD-based calculations and provide a 
much firmer basis compared to the quark model calculations for systematic evaluation of the uncer- 



tainties. The calculations used in the current results have been summarized nicely in Ref. 11. The 
LCSR are expected to be valid in the region ^ 16 GeV^. The light cone sum rules calculations 
use quark-hadron duality to estimate some spectral densities, and offer a "canonical" contribution 
to the related uncertainty of 10% with no known means of rigorously limiting that uncertainty. 
The theory community is currently d ebf^ ting the size of potential contributions to the form factors 
missing from the LCSR approachE2JlH3 that have been revealed using the newly-developed soft 
collinear effective theory (SCET). The B piv form factors, in particular, could be appreciably 
overestimated, biasing |Vuf,| low. Two exclusive results will therefore be presented in this review, 
one based on the full set of exclusive results, and the second based only on results in the > 16 
GeV^ region for B — > piu. 

The LQCD calculations that can be applied to experimental B — > X^ii^ decay remain, to 
date, in the "quenched" approximation (no light quark loops in the propagators), which limits the 
ultimate precision to the 15% to 20% range. The range accessible to these calculations has been 
<i 16 GeV^ Significant progress has been made towards unquenched lattice QCD calculations, and 



a recent comparisonllij of a range unquenched results to experiment shows much better agreement 
(few percent) than the corresponding quenched results. Work has begun on the unquenched form 
factors needed for \Vub\, though the initial results havp hppri limited to valence quarks closer to 



the strange quark mass. Nevertheless, initial results l^^^r^^ are compatible with the 15% to 20% 
uncertainties used for the quenching uncertainty, lending them some validity. 

The exclusive | Vub \ results are summarized in p?able 2.1 . Ttese include a simultaneous measure- 



ment of the B pi_7r£z^ and the_S piu transitions by CLEOlii^, and measurement of the B 



rate by CLEOtl^ and BaBaiti3. All measurements employ the missing energy and momentum 
to estimate the neutrino momentum. With that technique, the major background results from 
b civ decays in events that cannot be properly reconstructed (for example, because of additional 
neutrinos in the event) and hence which overestimate the neutrino energy. All measurements also 
employ the isospin relations 

r (S° ^ TT-tu) = 2T {B+ ^ vr"^ zy) 

and 

r {B" p-tv) = 2T {B+ p"e+u) 
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to combine the charged and neutral decays. These relationships can be distorted hy p — u> mixing 
and all results discussed here allow for this possibility in their systematic evaluation. 

In the combined vr and p measurement, strict event quality requirements were made that 
resulted in a low efficiency, but a relatively low background to signal ratio over a fairly broad lepton 
momentum range. The p-only analyses employ looser event cleanliness requirements, resulting in 
a much higher efficiency. The efficiency gain comes at the price of an increased background, and 
the analyses are primarily sensitive to signal with lepton momenta above 2.3 GeV/c, which is near 
(and beyond) the kinematic endpoint for b ciu decays which are therefore highly suppressed. 

The combined vr and p analysis of CLEO employs relatively loose lepton selection criteria and 
extracts rates independently in three separate intervals. Form factor dependence of the rates 
is then evaluated using models and calculations that exhibit a broad variation in dT/dq^, which 
shows that this approach has eliminated model dependence of the rates in niv, and significantly 
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Table 2.1: Summary of all exclusive \Vub\ measurements. For the CLEO '00 
and BaBar '01 measurements, the errors arise from statistical, experimental sys- 
tematic and form factor modeling uncertainties, respectively. For the CLEO '03 
measurements, the errors arise from statistical, experimental systematic, plu form 
factor, and LQCD and LCSR calculation uncertainties, respectively. In the CLEO 
'03 averages, the LQCD and LCSR uncertainties have been treated as correlated. 







mode 


IK.I(io-3) 




range 


FF 


CLEO 


'00 


piu 


3.23 zb 0.24l[!:l ± 0.58 




all 


model survey 


BaBar 


'01 


plv 


3.64 ± 0.22 ±0.25tg:^^ 




all 


model survey 


CLEO 


'03 


nil' 


3.33 ±0.24 ±0.15 ±0.06 




< 16 GeV^ 


LCSR 


CLEO 


'03 


nil' 


2.88 ± 0.55 ± 0.30 ± 0.18 ig j^ 




> 16 GeV2 


LQCD 


CLEO 


'03 


niv 


3.24 ± 0.22 ± 0.13 ± 0.09tEJ:i 




average 




CLEO 


'03 


piv 


2.67 ± 0.27 ± 0.17 t^Z 


q' 


< 16 GeV2 


LCSR 


CLEO 


'03 


piv 


3.34 ± 0.32 ± 0.47 


q' 


> 16 GeV^ 


LQCD 


CLEO 


'03 


piv 


3.00 ± 0.2ltl f, ± 0.28t°:^^ 




average 




CLEO 


'03 


n + p 


3.17 ±0.17tl]:i? ±0.0311]:^^ 




average 




CLEO 


'03 


n + p 


3.26 ± 0.19 ± 0.15 ±0.04t°j^ 




average 


no piu LCSR 



reduced it in piu. To further reduce modeling uncertainties, CLEO then extracts |I4b| using only 
the LQCD and LCSR QCD-based calculations restricted to their respective valid q'^ ranges, thereby 
eliminating modeling used for extrapolation. Averages of the CLEO results, with and without the 



low q region for piu are listed in [Table 2.1 . 



A more complete review of recent B —i- X^iv branching fractions—including analyses too 



incomplete for inclusion in this | Vub \ summaxji,, can be found in Referencdiisl. Of note is the recent 
evidence presented for B uoiv by BELLE 
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With all results resting on use of detector hermetiiity, the potential for significant correlation 
among the dominant experimental systematics existsiia. Results from the three measuremBiits 
have been averaged here assuming full correlation in these systematics. The piv-orAy resultsii^^tiS , 
which depend more heavily on modeling even for the LCSR and LQCD calculations, are deweighted 
by 5% in the average. This yields 

\Vub\ = (3.27 ± 0.13 ±0.19^°^^) X 10"^ 

where the errors arise from statistical, experimental systematic and form factor uncertainties, re- 
spectively. While similar in precision to the exclusive result in the previous \Vub\ minireview, this 
result relies much less heavily on modeling. Should the LCSR form factors prove to be overesti- 
mated, we also provide an average excluding any result using information for q^ < 16 GeV^ in the 
pii' modes, with the result 

iKbl = (3.26 ± 0.19 ± 0.15 ± 0.04+;;j^) x 10'^ 

where the errors arise from statistical, experimental systematic piv form factor uncertainties, and 
LQCD and LCSR (treated as correlated), respectively. 

The future for exclusive determinations of \Vub\ appears promising. Unquenched lattice calcu- 
lations are appearing, with very encouraging results. These calculations will eliminate the primary 
source of uncontrolled uncertainty in these calculations, and have already provided some validity 
to the quenching uncertainty estimate used in the results presented here. Simultaneously, the B 
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factories are performing very well, and very large samples of events in which one B meson has 
been fully reconstructed are already being used. This will allow a more robust determination of 
the neutrino momentum, and should allow a significant reduction of backgrounds and experimen- 
tal systematic uncertainties. The high statistics should also allow more detailed measurements of 
dT/dq^, which have already provided a sorely-needed litmus test for the form factor calculations 
and reduced the form factor shape contribution to ±he uncertainty on | Vub \ ■ Should theory allow 



use of the full range of in the extraction of |14b|lH3, the B factories have already logged data 
sufficient for a 5% statistical determination of \Vub\- 

For both lattice and the B factories, nil' appears to be a golden mode for future precise 
determination of \Vub\- The one caveat is management of contributions from the B* pole, but recent 
workCa suggests that this problem can be successfully overcome. B — > rjiu will provide a valuable 
cross-check. The piv mode will be more problematic for high precision: the broad width introduces 
both experimental and theoretical difficulties. Experiments must, for example, assess potential 
nonresonant vrvr contributions, but only crude arguments based on isospin and quark-popping have 
been brought to bear to date. Theoretically, no calculation, including lattice, has dealt with the 
width of the p. When the lattice calculations become unquenched, the p will become unstable 
and the vrvr final state must be faced by the calculations. The methodology for accommodation of 
high-energy two particle final states on the lattice has yet to be developed. The ujiiy mode may 
provide a more tractable alternative to the p mode because of the relative narrowness of the a; 
resonance. Agreement between accurate \Vub\ determinations from Triu and from Loii^ will provide 
added confidence in both. 



3. Combined results 

The experimental bounds provided for the outstanding uncertainties in the inclusive \Vub\ mea- 
surements in the low Mx, high region and theoretical work which clarifies the reliability of 
the LCSR and quenched LQCD form factors make possible comparison of the inclusive and exclu- 
sive determinations of \Vub\- Results agree to better than 1.5 times the quadratic combination of 
the quoted uncertainties. Therefore it becomes feasible to propose an average the inclusive and 
exclusive results, which have comparable accuracies. 

The uncertainties have been combined in quadrature, using the larger (upward) error for the 
exclusive numbers. The proposed average of the inclusive and exclusive results, with all the exclusive 
data considered, is 

iKbl = (3.67 ±0.47) X 10"^ 

Including in the average only the exclusive analyses based on data with > 16 GeV^ in the piv 
mode, the average becomes \Vub\ = (3.70 ± 0.49) x 10~^, so exclusion of this region, if appropriate, 
has only a minor effect. 

The procedure proposed here results in a value of \Vub\ with a 13% uncertainty. With the 
experimental and theoretical progress expected over the next few years an improvement of the 
accuracy at the 10% level, and possibly below, appears now realizable. 
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